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It has recently been pointed out that the 511 keV emission line detected by 
INTEGRAL/SPI from the bulge of our galaxy could be explained by annihilations 
of light dark matter particles into e + e~ . We present the possibility that dark matter 
could be made of scalar candidates, namely , of the Higgs bosons in the models based 
on SU(3)c <8> SU(3)l (g> U(l)jv (3-3-1) gauge group. These particles are singlet of the 
SU(2)l (g) U(l)y group, so they do not interact with the ordinary particles, exept 
the Higgs boson in the standard model. The Spergel-Steinhardt condition for self- 
interacting dark matter gives a bound on the mass of the candidates to be a few 
MeVs. Besides the scalar candidates, which exist in both non-SUSY and SUSY 
3-3-1 models with right-handed neutrinos, the spin ^ candidate exists in a variant 
3-3-1 version with exotic neutral lepton. In contrast to the singlet models, where an 
extra symmetry must be imposed to account the stability of the dark matter, here 
the decay of the candidates is automatically forbidden in all orders of perturbative 
expansion. This is because of the following feature: these scalars are singlets, i.e., 
in bottom of the Higgs triplet. Therefore, the standard model fermions and the 
standard gauge bosons cannot couple with them. 



1 Introduction 

Observation of the cosmic microwave background, the primordial abundances of light 
elements and large scale structure have revealed that a great deal of the mass of our 
universe consists of dark matter (DM) . It is an amazing fact that even as our under- 
standing of cosmology progresses by leaps and bounds, we remain almost completely 
ignorant about the nature of most of the matter in the universe [1]. It has recently 
been pointed out that the 511 keV emission line detected by INTEGRAL/SPI from 
the bulge of our galaxy could be explained by annihilations of light dark matter 
particles into e + e~ . The nature of DM is still a challenging question in cosmology. 
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Cosmological models with a mixture of roughly 35% collisionless cold DM such as 
axions, WIMPs, or any other candidate interacting through the weak and gravi- 
tational forces only, and 65% vacuum energy or quintessence match observation of 
the cosmic microwave background and large scale structure on extra-galactic scales 
with remarkable accuracy [2,3]. It is known that only a fraction of the dark matter 
can be made of ordinary baryons and its enormous amount has unknown, nonbary- 
onic origin [4]. Until a few years ago, the more satisfactory cosmological scenarios 
were those ones composed of ordinary matter, cold DM and a contribution associ- 
ated with the cosmological constant. To be consistent with inflationary cosmology, 
the spectrum of density fluctuations would be nearly scale-invariant and adiabatic. 
However, in recent years it has been pointed out that the conventional models of 
collisionless cold DM lead to problems with regard to galactic structures. They 
were only able to fit the observations on large scales (S> 1 Mpc). Also, iV-body 
simulations in these models result in a central singularity of the galactic halos [5] 
with a large number of sub-halos [6], which are in conflict with astronomical ob- 
servations. A number of other inconsistencies are discussed in Refs. [7,8]. Thus, 
the cold DM model is not able to explain observations on scales smaller than a 
few Mpc. It has recently been shown that an elegant way to avoid these prob- 
lems is to assume the so called self-interacting dark matter [9]. One should notice 
that, in spite of all, self-interacting models lead to spherical halo centers in clusters, 
which is not in agreement with ellipsoidal centers indicated by strong gravitational 
lens observations [10] and by Chandra observations [11]. However, self-interacting 
dark matter models are self- motivated as alternative models. It is a well-accepted 
fact that the plausible candidates for DM are elementary particles. The key prop- 
erty of these particles is that, they must have a large scattering cross-section and 
negligible annihilation or dissipation. The Spergel-Steinhard model has motivated 
many follow-up studies [4, 12, 13] . Several authors have proposed models in which 
a specific scalar singlet that satisfies the self-interacting dark matter properties is 
introduced in the standard model (SM) in an ad hoc way [4,13]. The SM offers no 
options for DM. The first gauge model for SIDM is the extended SM based on the 
SU(3) c ®SU(3) l ®U(1)at (3-3-1) gauge group [14,15]. The 3-3-1 model founded by 
Fregolente and Tonasse [16] to be SIDM differs from another one founded by Long 
and Lan [17]. Models with SU (3)c <8> SU '(3) l <8> U (1) n gauge symmetry (called 3-3-1 
models for short) are interesting possibilities for the physics at the TeV scale [14,15]. 
At low energies they coincide with the standard model and some of them give at 
least partial explanation to some fundamental questions that are accommodated but 
not explained by the standard model. For instance, in order to cancel the triangle 
anomalies, together with asymptotic freedom in QCD, the model predicts that the 
number of generations must be three and only three; (ii) the model of Ref. [14] 
predicts that (g'/g) 2 = s'm 2 9\v/(l — 4sin 2 #w0> thus there is a Landau pole at the 
energy scale \i at which s'm 2 6yy(^) = 1/4. According to recent calculations fj, ~ 4 
TeV [18, 19] ; (Hi) the quantization of the electric charge [20] and the vectorial 



Candidates for Dark Matter 



3 



character of the electromagnetic interactions [21] do not depend on the nature of 
the neutrinos i.e., if they are Dirac or Majorana particles; (iv) as a consequence of 
item ii) above, the model possesses perturbative M = 1 supersymmetry naturally 
at the fj, scale [22,23]; (v) the Peccei-Quinn [24] symmetry occurs naturally in these 
models [25] ; (vi) since one generation of quarks is treated differently from the others 
this may be lead to a natural explanation for the large mass of the top quarks [26]. 
Moreover, if right-handed neutrinos are considered transforming non-trivially [15], 
3-3-1 models [14] can be embedded in a model with 3-4-1 gauge symmetry in which 
leptons transform as (ui,l,uf, I c )l ~ (1,4,0) under each gauge factors [27,28]. The 
SU(3)l symmetry is possibly the largest symmetry involving the known leptons 
(and SU(4:)l if right-handed neutrinos do really exist). This make 3-3-1 or 3-4-1 
models interesting by their own, and it has been the source of interest recently [29] 
because it requires that the number of fermion families be a multiple of the quark 
color in order to cancel anomalies, which suggest a path to the solution of the flavor 
problem. The Peccei-Quinn [24] symmetry occurs naturally in these models [25]. 
Another important feature of these models is that the SU(2)^ group is totally em- 
bedded in SU(3)l. A subject that has not been given much attention by particle 
physicists in the past, could prove to be a remarkable powerful and precise probe 
of the properties of dark matter. This paper is organized as follows: In Sec. 2 we 
briefly introduce necessary elements of the non-SUSY 3-3-1 models, and the can- 
didates for dark matter. Sec. 3 is devoted to supersymmetric 3-3-1 model with 
right-handed neutrinos. We summarize our result and make conclusions in the last 
section. The aim of this reviewr is to show that the 3-3-1 models contains can- 
didates for dark matter. The outline of this paper is as follows. In Sec. 3, the 3-3-1 
models are briefly recalled. Properties of DM are introduced. The candidates for 
SIDM in four 3-3-1 models are explored. The last section is devoted for conclusions. 



To frame the context, it is appropriate to recall briefly some relevant features of the 
3-3-1 models [14,15]. We first introduce the minimal version proposed by Pisano, 
Pleitez and Frampton [14] 

2.1 The 3-3-1 model with exotic lepton 

The model treats the leptons as SU(3) L triplets [30] 



where o = 1,2,3 is the generation index. Each charged left-handed fermion field 
has its right-handed countepart transforming as a singlet of the SU(3)l group. Two 



2 The non-SUSY 3-3-1 models 
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of the three quark generations transform as antitriplets and the third generation is 
treated differently - it belongs to a triplet: 



QiL 




(3,3, 



1, 

3' 



(2) 



u iR ~ (3, 1, 2/3), d lR ~ (3, 1, -1/3), Afl ~ (3, 1, -4/3), i = 1, 2, 

Qsl= d 3 L ~ (3,3,2/3), (3) 

V ^ / 

u 3 R ~ (3, 1, 2/3), d 3 i? ~ (3, 1, -1/3), T B ~ (3, 1, 5/3). 

Of the nine gauge bosons W a {a = 1,2, 8) and B of SU(3)l and U(1)n, four are 
light: the photon (A), Z and VF ± . The remaining five correspond to new heavy 
gauge bosons Zi, and the doubly charged bileptons X ±=t . Symmetry breaking 
and fermion mass generation can be achieved by three scalar SU(3) L triplets 

SU(3) C ® SU(3) l ®U(1)jv 

i(x) 

SU(3) C ® SU(2) L ®U(l)y 

i<»7>,<P> 
SU(3)c7 ® U(1) Q , 

where the minimally required scalar multiplets are summarized as 

[ X ~ 

X = X 




(4) 

The vacuum expectation value (VEV) (x T ) = (0,0, w/V^) yields masses for the 
exotic quarks, the heavy neutral gauge boson (Z2) and two charged gauge bosons 
(X ++ ,Y + ). The masses of the standard gauge bosons and the ordinary fermions 
are related to the VEVs of the other scalar fields, (r)°) = v/y/2, (p°) = u/y/2. In 
order to be consistent with the low energy phenomenology we have to assume that 
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u v, u). By matching the gauge coupling constants we get a relation between g 
and <7at - the couplings associated with SU(3) L and U(l)jy, respectively: 

g% 6 s 2 w {M Z2 ) 

? l-4s* w (M Z2 y W 

where e = g sw is the same as in the SM. The most economical, gauge invariant 
and renormalizable Higgs potential is 

V(r),P,X) = Pirfv + f4p*P + t4x*X + a i (v*v) + a 2 (ftp) +a 3 (x < x) 
(V ??) (aAP ] p + a 5 x^x) + a 6 (p 1 p) (x'x) + a 7 (p^v) (Vp) 
-as (x j v) (V x) + a 9 (p f x) (x ] p) + \ (f^ViPjXk + H. c.) .(6) 



+ 
+c 



Here the p's and / are coupling constants with dimension of mass with 03 < and 
/ < from the positivity of the scalar masses [31]. 

Symmetry breaking is initiated when the scalar neutral fields are shifted as 
f = Vtp + £tp + i(tp, with tp = rj°, p°, x°- The details of the physical spectrum of 
the neutral scalar sector are crucial for our results. It is given in Refs. [31] and we 
summarize them here. Firstly we notice that real part of the shifted fields leads to 
the three massive physical scalar fields Hf, H®, defined by 



« ^ (7) 



where we are using w 3> v,u. The scalar H® is the one that we can identify with 
the standard model Higgs, since its squared mass, 

2 _ , a2U 4 - aiv A 
m i ~ 4 2 a - ' V 8 

carries no any feature from the 3-3-1 breakdown to the standard model. On the 
other hand H®, with squared mass 

ml » -4a 3 w 2 , (9) 

is a typical 3-3-1 scalar. So, there is no any massless Goldstone boson rising from 
the real part of the neutral sector. On the other hand, from the imaginary part we 
have two Goldstone and one massive physical state h° with eigenstate 

Cx « h° (10) 

and squared mass 



vuw 
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It is important to notice that and Q p are pure massless Goldstone states. The 
approximation in Eqs. ©, ©, © and (fTU|) . is valid for w S> v,u. This condition 
leads to relations among the parameters of the scalar potential ([U]). One of them, 
which enters in the H^hPhP interaction, is 

a 5 v 2 + 2a 6 u 2 « — — (12) 

We must consider also the matter coupling through the scalar fields. In the model 
of Ref. [30], the full Yukawa Lagrangians that must be considered are 

C t = "Ef^^^ 

ab V/ 7 



£q = QiLT,( G{ i J b )u mV + G^ ) D bR p + G^J mX ) + 



E QotL f?W + + E Wj + H - c - ( 14 ) 

where C$, Gf h \ C$\ F%\ and F%j are the Yukawa con- 

pling constants, rf , p* and x* denote the rj, p and X antiparticle fields, respectively. 
The main properties that a good dark matter candidate must satisfy are stability 
and neutrality. Therefore, we go to the scalar sector of the model, more specifically 
to the neutral scalars, and we examine whether any of them can be stable and in 
addition whether they can satisfy the self- interacting dark matter criterions [9]. In 
addition, one should notice that such dark matter particle must not overpopulate 
the Universe. On the other hand, since our dark matter particle is not imposed 
arbitrarily to solve this specific problem, we must check that the necessary values 
of the parameters do not spoil the other bounds of the model. 

We can check through a direct calculation by employing the Lagrangians ([6]), 
(|13p and (|14p and by using the eigenstates ([7]) and (|10p that the Higgs scalar hP and 

can, in principle, satisfy the criterions above. Remarkably they do not interact 
directly with any standard model field except for the standard Higgs H®. However, 
h° must be favored, since we have checked that it is easier to obtain a large scattering 
cross section for it, relative to by a convenient choice of the parameters. In 
contrast to the singlet models of the Refs. [4, 13], where an extra symmetry must 
be imposed to account the stability of the dark matter, here the decay of the hP 
scalar is automatically forbidden in all orders of perturbative expansion. This is 
because of the following features: i) this scalar comes from the triplet %, the one 
that induces the spontaneous symmetry breaking of the 3-3-1 model to the standard 
model. Therefore, the standard model fermions and the standard gauge bosons 
cannot couple with h°, ii) the h scalar comes from the imaginary part of the Higgs 
triplet x- As we mentioned above, the imaginary parts of r/ and p are pure massless 
Goldstone bosons. Therefore, there is not physical scalar fields which can mix with 
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h°. So, the only interactions of h° come from the scalar potential and they are 
H$h°h° and H^h°h°. This latter has strength 2i (a 5 v 2 + a 6 u 2 ) /v w = 2iQ. We can 
check also that h° does not interact with other exotic particle. 

Hence, if v ~ u ~ (100 - 200) GeV and -1 < a 5 ~ a 6 < 1, the h° can 
interact only weakly with ordinary matter through the Higgs boson of the stan- 
dard model H®. The relevant quartic interaction for scattering is h°h h°h , whose 
strength is —ia^. Other quartic interactions evolving h° and other neutral scalars 
are proportional to 1/w and so we neglect them. The cross section of the pro- 
cess 

h o h o ^ h o h o 

via the quartic interaction is a = a^/GAiTm^. The contribution 
of the trilinear interactions via H® and H® exchange are negligible. There is no 
other contribution to the process involving the exchange of vector or scalar bosons. 
A self- interacting dark matter candidate must have mean free path A = 1/na in 
the range 1 kpc < A < 1 Mpc, where n = p/m^ is the number density of the 
h° scalar and p is its density at the solar radius [6]. Therefore, with 03! = — 1, 
-0.208 x 10~ 7 GeV < / < -0.112 x 10~ 6 GeV, w = 1000 GeV, u = 195 GeV 
and p = 0.4 GeV/cm 3 , we obtain the required Spergel-Steinhardt bound, i. e., 
2 x 10 3 GeV~ 3 < a/m h < 3 x 10 4 GeV^ 3 [9]. 

With this set of parameter values, we see from Eq. (Ilip that 5.5 MeV < < 29 
MeV. This means that our dark matter particle is non-relativistic in the decoupling 
era (decoupling temperature ~ 1 eV) and, for a standard model Higgs boson mass 
~ 100 GeV [32], it is produced by a thermal equilibrium density of the standard 
Higgs scalar to h°h° pairs [13]. The density of the h° scalar from the decay can 
be obtained following the standard procedure, i. e., we must solve the Boltzmann 
equation 

^± + 3Hn h = (T H )np\ (15) 

where nh is the number density of the h° scalar at the time t, H is the Hubble 
expansion rate, 

r » = si (16> 



is the decay rate for the with energy E and 



,(eq) 



1 roo EJE 2 



mr 



2vr2 J m , e E / T - 1 



(17) 



is the thermal equilibrium density of the standard Hf at temperature T [33]. We 
are using the condition that the temperature is less than the electroweak phase 
transition Tew > 1.5mi [13]. The thermal average of the decay rate is given by 



(18) 
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(19) 



where a is an integration parameter that can be taken to be 1.87 [13]. We define 
/3 = nh/T 3 and in the radiation-dominated era we write the evolution equation (|15p 

as 

dT KT 3 87r 3 Ke m ^ T \T 2 J 

where K 2 = 4-7r 3 g (T) /45mpj, (3^ cq -^ = n h cq ^ /T 3 is the (3 parameter in the thermal 
equilibrium, mp\ = 1.2 x 10 19 GeV is the Planck mass and g (T) = qb + 7o,f/8 = 
136.25 for the model of the Ref. [30]. gs and gF are the relativistic bosonic and 
fermionic degrees of freedom, respectively. Here we are taking T = mi since this 
regime gives the larger contribution to /3 [13]. Hence, 

aG 2 

" = S«5?- <20> 

Now, the cosmic density of the h° scalar is 

n h = 2g(T r )T 3 ^-, (21) 



where T 7 = 2.4 x 10 -4 eV is the present photon temperature, g (T 7 ) = 2 is the 
photon degree of freedom and p c = 7.5 x 10 -47 /i 2 , with h = 0.71, being the critical 
density of the Universe. Let us take mh = 7.75 MeV, v = 174 GeV, as = 0.65, 
— ciq = 0.38 (actually in our calculations, we have used a better precision for as 
and o,q) and mi = 150 GeV. Thus, from Eqs. ([20]) and (f2T1) we obtain Qh = 0.3. 
Therefore, without imposing any new fields or symmetries, the 3-3-1 model possesses 
a scalar field that can satisfy all the properties required for the self-interacting dark 
matter and that does not overpopulate the Universe. 

The candidate for self-interacting dark matter that we propose here differs from 
the singlet models of Refs. [4, 13] in an important point. As we have discussed 
above it comes from a gauge model proposed with another motivation that has an 
independent phenomenology. Therefore, the values of the parameters that we impose 
here must not spoil the preexisting bounds. We can obtain mi ~ 150 GeV from 
Eq. ([8]) with ai = 1.2 and a2 = 0.36. From Eq. ([9]) we have m% = 1 TeV. On the 
other hand, one should notice that mh has a small value since — / ~ (10~ 7 — 10 -6 ) 
GeV and u ~ 195 GeV. However, hP does not couple to the particles of the standard 
model for the Higgs boson. Thus, it evades the present accelerator limits. The 
constants as and do not enter in the masses of the particles of the model and so, 
it is free in this work [31]. 

2.2 The 3-3-1 model with right-handed neutrinos 

In this model the leptons are in triplets, and the third member is a RH neutrino: 

ft = K, eh KD T ~ (1, 3, -1/3), 4 ~ (1, 1, -1). (22) 
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Each charged left-handed fermion field has its right-handed counterpart transform- 
ing as a singlet of the SU(3)l group. The first two generations of quarks are in 
antitriplets while the third one is in a triplet: 

Q l L = (d iL ,-u iL ,D iL ) T ~ (3,3,0), (23) 

u iR ~ (3, 1, 2/3), d lR ~ (3, 1, -1/3), Afl ~ (3, 1, -1/3), i = 1, 2, 

Qsl = («ai, dsi, T L ) T ~ (3, 3, 1/3), (24) 
u 3fi ~ (3, 1, 2/3), d 3R ~ (3, 1, -1/3), Tr ~ (3, 1, 2/3). 
The charged gauge bosons are defined as 

V2W+ = Wl-iW*,V2Y- = W*-iW*, 
V2X° = Wj-iW* (25) 

The physical neutral gauge bosons are again related to Z, Z' through the mixing 
angle <f>. The symmetry breaking can be achieved with just three SU(3) l triplets 



X = 


(x° 


<x 


,X'°) T ~ (1,3,-1/3), 




p = 




,P°; 


P' + ) T ~ (1,3,2/3), 


(26) 


V = 






l 7 ? '°) T ~ (1,3,-1/3),. 





The necessary VEVs are 

(X) = (0, 0, u/V2) T , (p) = (0, u/y/2, 0) r , (r?) = (v/y/2, 0, 0) T . (27) 
After symmetry breaking the gauge bosons gain masses 

m 2 w = \g 2 (u 2 + v 2 ), M 2 =\ 9 2 (v 2 + u 2 ),M 2 = 1 - g 2 (u 2 + u> 2 ). (28) 

Eqn. (f58l) gives us a relation 

v ^ = u 2 + v 2 = 246 2 GeV 2 . (29) 

In order to be consistent with the low energy phenomenology we have to as- 
sume that (x) 3> (p), (rj) such that mw ^ Mx,My- The symmetry-breaking 
hierarchy gives us splitting on the bilepton masses [34] 

\M X -M Y \ < m 2 w . (30) 

Our aim in this paper is to show that the 3-3-1 model with RH neutrinos furnishes 
a good candidate for (self-interacting) dark matter. The main properties that a 
good dark matter candidate must satisfy are stability and neutrality. Therefore, we 
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go to the scalar sector of the model, more specifically to the neutral scalars, and 
we examine whether any of them can be stable and in addition whether they can 
satisfy the self-interacting dark matter criterions [9]. In addition, one should notice 
that such dark matter particle must not overpopulate the Universe. On the other 
hand, since our dark matter particle is not imposed arbitrarily to solve this specific 
problem, we must check that the necessary values of the parameters do not spoil the 
other bounds of the model. Under assumption of the discrete symmetry x — X> 
the most general potential can then be written in the following form [35] 

V(v,P,X) = PiV + V + P2P + P + pIx + X + Ai(r? + r?) 2 + X 2 (p + p) 2 + A 3 (x + x) 2 

+ (v + v)[Mp + p) + A 5 (x + x)] + a 6 (p + p)(x + x) + \ 7 (p + v)(v + p) 

+ A 8 (x + r?)(r/ + x) + A 9 (p + x)(x + p) + Aio(x + ^ + ?? + x) 2 - (31) 
We rewrite the expansion of the scalar fields which acquire a VEV: 

rf = -j= (v + 6? + K v ) ; p° = (« + t P + K P ) ; x° = ^= (w + £ x + K x ) • (32) 

For the prime neutral fields which do not have VEV, we get analogously: 

v° = ^ (?; + <); x'° = 7^ (& + <)■ (33) 

Requiring that in the shifted potential V, the linear terms in fields must be absent, 
we get in the tree level approximation, the following constraint equations: 

\ i 

p\ + X x v 2 + -X 4 u 2 + -X 5 w 2 = 0, 

p 2 2 + X 2 u 2 + ^X A v 2 + ^X 6 w 2 = 0, (34) 

fij + X 3 w 2 + -X 5 v 2 + -X 6 u 2 = 0. 

Since dark matter has to be neutral, then we consider only neutral Higgs sector. 
In the ^rj, £p, £ x , , basis the square mass matrix, after imposing of the constraints 
(|34p . has a quasi-diagonal form as follows: 



where 



m! h - - 



j / 2X±v 2 X4VU X5VW 

X4VU 2X2U 2 Xquw I , (36) 
* X§vw Xquw 2A3U) 2 



and 

M 2 H =[ C!7) 
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The above mass matrix shows that the prime fields mix themselves but do not mix 
with others. In the limit 

X±v , X2U, A4U <C X5W, Xqw, (38) 
we obtained physical eigenstates H\(x) and cr(x) 

( Hx(x) \ 1 ( A 6 n -X 5 v \( Z v \ 



\ a(x) J (X 2 5 v 2 + X 2 6 u 2 ) 1 / 2 \ X 5 v X 6 u J \ £ p 
with masses [35] 



(39) 



2 2 
m 2 Hl ~ ^ r (2AiA 6 -A 4 A 5 ) « -^-(2A 2 A 5 - A 4 A 6 ), (40) 

2 _ 1, 2 . A 4 A 6 n 2 1 2 A4A5U 2 
m CT « -An, + « -A 2 . + (41) 

Eqs. (|40p and (|4ip also give us relations among coupling constants and VEVs. An- 
other massive physical state H$ with mass: 



■2 



13 



-A 3 ^ 2 . (42) 



The scalar cr(x) is the one that we can identify with the SM Higgs boson [35]. In 
the approximation w 3> v, mass matrix M 2 H gives us one Goldstone £^ and one 
physical massive field with mass 

m| = -(^ + A 10 )^. (43) 

In the pseudoscalar sector, we have three Goldstone bosons which can be identified 
as follows: G2 = C?7> G3 = ( p , G4 = C, x and in the ('°, ('° basis 

A« \ I w 2 vw 



We easily get one Goldstone G' 5 and one massive pseudoscalar boson (L with mass 



m?,=-(^ + A 10 )«, 2 . (45) 

It is to be emphasized that, both £' and Q' are in an singlet of the SU (2). Therefore 
they do not interact with the SM gauge bosons W^,Z° and 7. Unlike the 3-3-1 
model considered in [16], here we have two fields which can be considered as dark 
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matter. To get the interaction of dark matter to the SM Higgs boson, we consider 
the following relevant parts 



+ \x, 



v 2 + 2v£ ri + e 2 + C' + Q + Q + 2?? V 



v 2 + 2v^ + £ 2 + Q 2 + Q + Q 2 + 277+77' 



x n 2 + 2< p + ^ + Cp + 2p + p-+2p-p 
Substituting (|39p we get couplings of SIDM with the SM Higgs boson o 



(46) 



u(x) 



X^ + X 2 



u- 



G 2 (X) 



2{x 2 y + xiv?) 



X x X b v" 



xiv + 



XdXf 



+ 



Hi(x)a{x) 



(A|u 2 + X 2 u 2 



Xi — -77 I X 5 X 6 uv 



t'2 + C '2 



(47) 



From Yukawa couplings, we see that our candidates do not interact with ordinary 
leptons and quarks [36]. 



^Yuk 



X3aQ3LUaRV + X iia Q iL d aR r]* + h.c. 

X 3a (u3Lri + dzLTf + f L 7]'°)u aR + A 4ia (d iL rf* - U iL T] + + D iL T)'°*)daR + h.c. 



We see that the candidates for dark matter in this model have not couplings with 
all the SM particles except for the Higgs boson. For stability of DM, we have to 
put mass of the SM Higgs boson is twice bigger mass of the candidate 



1 2 A 4 A 6 u 2 
2 XlV + ^AT 



1 2 A4A5W 2 

2 X2U + ^AT " 



+ A10 )w 



To avoid the interaction of DM with Goldstone boson, we have 

A - A4 

Al "T 



(48) 



(49) 



The wrong muon decay (/i — > e v e v^) gives a lower limit for singly charged 
bilepton My ~ 230 GeV. Combining Eqns. Mfy with ([60} we obtain the 

following relation: u ~ v w 100 - 200 GeV and w « (500 - 1000) GeV. The 
cross section for hh — > hh (where h stands for £' and with quartic interaction is 
a = X\/A-Km 2 l . The requirement on the quality o~ e i/(mh[GeV\) denoting the ration 
of the DM elastic cross section to its mass (measured in GeV) is that [9, 13,37] 



2.05 x 10 3 GeV" 3 < < 2.57 x 10 4 GeV" 3 



(50) 
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Taking Ai = 1 we get 4.7 MeV < m h < 23 MeV. The SIDM candidates interact 
with the SM Higgs boson by strength 0.65 if A5 = Ag = 1 and u = v = 175 GeV are 
taken. Now consider the cosmic density of the h scalar given by [16]: 

- 2 ' AT ^^m- (51) 

where T 7 = 2.4 x 10 -4 eV is the present photon temperature, g(Tj) = 2 is the 
photon degree of freedom and p c = 7.5 x 10~ 47 /i 2 with h = 0.71, being the critical 
density of the Universe. Taking = 4.7 MeV, we obtain = 0.18. This means 
that the SIDM candidates do not overpopulate the Universe. 



2.3 The 3-3-1 model with exotic neutral lepton 

In this model the leptons are in triplets, and the third member is a neutral exotic 
lepton and left-handed fermion field has its right-handed counterpart transforming 
as a singlet of the SU(3)l group [38]: 

ft = K, e a L , Ntf ~ (1, 3, -1/3), e a R ~ (1, 1, -1), JVg ~ (1, 1,0). (52) 

The first two generations of quarks are in antitriplets while the third one is in a 
triplet: 

QiL = (d iL ,-u iL ,D iL ) ~ (3,3,0), (53) 
u iR ~ (3, 1, 2/3), d lR ~ (3, 1, -1/3), Afl ~ (3, 1, -1/3), i = 1, 2, 

Qzl = (u 3L , d 3L ,T L ) T ~ (3, 3, 1/3), (54) 
u 3R ~ (3, 1, 2/3), d 3R ~ (3, 1, -1/3), T R ~ (3, 1, 2/3). 
The charged gauge bosons are defined as 

\/2W+ = Wl-iWl,y/2Y~ = Wl-iWl t 

= Wi-iWl. (55) 

The physical neutral gauge bosons are again related to Z, Z 1 through the mixing 
angle <fi. The symmetry breaking can be achieved with three SU(3)l triplets and an 
sextet 

X = (x°,X",X'°) T ~ (1,3,-1/3), 

p = (p+,p°,p> + ) T ~ (1,3,2/3), (56) 



Candidates for Dark Matter 



14 



The necessary VEVs are 

( X ) = (0, 0, lo/V2) t , (p) = (0, u/y/2, 0) r , (n) = (v/y/2, 0, 0) T . (57) 
The VEV of sextet has the form 





( 





Va 2 


(S) = 













V V «2 









After symmetry breaking the gauge bosons gain masses 

m 2 w = ^g 2 (u 2 +v 2 +v A ), M 2 =-^g 2 (v 2 +uj 2 +v 2 ),M 2 x = - A g 2 (u 2 +u J 2 +v 2 ). (58) 

Eqn.(|5"8"]) gives us a relation 

v ^ = u 2 + v 2 + v \ = 246 2 GeV 2 . (59) 

In order to be consistent with the low energy phenomenology we have to as- 
sume that (x) S> (p), (77) such that mvi/ <C Mx,My- The symmetry-breaking 
hierarchy gives us splitting on the bilepton masses [34] 

\M 2 X -M§\ < m 2 w . (60) 

The Yukawa couplings are given 

C Y = fiL a 4>u R + ^fiL a SL b (61) 

The mass eigenstates (u e ,N) are related to the weak eigenstates (u e ,N) 

v = v cos 9 - N sin 9, N = v sin 9 + N cos 9 (62) 

where 9 is a mixing angle. A mass distribution extended well beyond the visible 
galaxy [39] 




Fig.l: Two diagrams for the decay N — > v e ^ 
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<? 7 

H ^ # 

* i * i „ 

v e e f e 

Fig. 2: Diagram for Higgs contribution to the decay N — > i/ e 7 

The decay of iV particle into neutrino and photon contains three graphs depicted in 
Fig.l which give [38] 

T ^ N -^=TM^ W{2e) (64) 

which corresponds to lifetime 

r « 4.67 x 10 14 sin" 2 (20)(lkeV/m7vc 2 ) 5 yr (65) 

It is obvious that the N- particle's life can be longer than the age of the universe, 
if rriM ~ 1 kEV. It was shown that the Higgs contribution is much smaller than 
those of the two first, and the N- particle meets the constraints on dark matter 
from cosmology and galaxy formation. 

3 Supersymmetric 3-3-1 model with right-handed neu- 
trinos 

Supersymmetric version on the 3-3-1 model with right-handed neutrinos has been 
proposed in [22,40]. Here we will follow the usual notation writing for a given 
fermion /, the respective sfermions by / i.e., I and q denote sleptons and squarks 
respectively. Then, we have the following additional representations 

fu 3 \ 

d 3 -(3,3,1/3), 

(66) 



QaL 



L a L 



( d a \ 

u a ~ (3,3*,0), Q 3L 

fi>a\ 

l a -(1,3,-1/3), 
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(1,1,1) 




(3*,1 



2/3), d\ L J c aL 



1,1/3) 



(67) 



with a = e, u,,t; i = 1,2,3; and a = 1,2. However, when considering quark (or 
squark) singlets of a given charge we will use the notation u c iL ,d c iL (v,iL,df L with 
= 1,2,3. The supersymmetric partner of the scalar Higgs fields, the higgsinos, 

are 



and the respective extra higgsinos, needed to cancel the chiral anomaly of the hig- 
gsinos in Eq. (f68j) . are 



and the corresponding scalar partners denoted by r/,x', p' , with the same charge as- 
signment as in Eq. and with the following VEVs: v' = )/v / 2, w' = (x'^/Vz 
and v! = ( / o /0 )/v / 2- This complete the representation content of this supersymmet- 
ric model. Concerning the gauge bosons and their superpartners, if we denote the 
gluons by g b the respective superparticles, the gluinos, are denoted by A^, with 
b = 1, . . . , 8; and in the electroweak sector we have V b , the gauge boson of SU (3)l, 
and their gauginos partners A^; finally we have the gauge boson of U(1)n, denoted 
by V', and its supersymmetric partner A#. This is the total number of fields in the 
minimal supersymmetric extension of the 3-3-1 model of Refs. [15,30]. 

3.1 Superfields 

The superfields formalism is useful in writing the Lagrangian which is manifestly 
invariant under the supersymmetric transformations [41] with fermions and scalars 
put in chiral superfields while the gauge bosons in vector superfields. As usual the 
superfield of a field 4> will be denoted by <\> [42]. The chiral superfield of a multiplet 
4> is denoted by 




(68) 




(69) 




Candidates for Dark Matter 



17 



+ V2 0<f>(x) + ^=90 Bo m d m 4>{x) 
a/2 

+ OdF <t> {x), (70) 
while the vector superfield is given by 

v( x , e, o) = -0a m ev m {x) + mBv(x) - mev{x) 

+ -66dOD{x). (71) 

The fields F and D are auxiliary fields which are needed to close the supersymmetric 
algebra and eventually will be eliminated using their motion equations. For fermion 
superfields we use the notation 

LaL, Vol, QaL, Q?,L, d c iL , u'£, d'£ L . (72) 

For scalar superfields we write: fj, x, p and similar expressions for fj' , x' , p' and we 
must change (field) by (field)'. The vector superfield for the gauge bosons of each 
factor SU(3)c, SU(3)l and J7(1)jv are denoted by Vc,V& V,V; and V, respec- 
tively where_we have defined % = T b V£, V = T b V b ; V c = f b V b c , V = f b V b ; 
T b = X b /2, T b = —X* b /2 are the generators of triplet and antitriplets representa- 
tions, respectively, and X b are the Gell-Mann matrices. The Lagrangian of the model 
has the following form 

£3315 = &SUSY + £ so ft' (73) 

where Csusy is the supersymmetric part and C so ft the soft terms breaking explicitly 
the supersymmetry. 

3.2 The supersymmetric Lagrangian 

The supersymmetric part of the Lagrangian is decomposed in the lepton, quark, 
gauge, and the scalar sectors as follow: 

Csusy = ^Lepton + £ Quark + £ Gauge + ^Scalar ( 74 ) 

where 

^Leptons = [ W 2 ^* ' U + , (75) 

in the lepton sector, we have omitted the sum over the three lepton family for 
simplicity, and 



^Quarks = ' d 9 



(76) 
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in the quark sector, and we have denoted g s ,g,g' the gauge coupling constants for 
the SU(3)c, SU(3)l,U(1)n factors, respectively. In the gauge sector we have 



^Gauge = \j d2 ° t W cWc + WW + WW } 

- J d 2 e [W C W C + ww + WW ] 



+ 



where Wc, W e W are fields that can be written as follows 



W cc = -J-DDe 



W ? 



8g s 
1 



8g 
1 



DDD(V\ ( = 1,2. 



Finally, in the scalar sector we have 



C 



dH 



Escalar 



+ / d z 9W + / d z 6W, 
where W is the superpotential 



2/3W i / j2/ 



(77) 



(78) 



(79) 



3.3 The scalar Potential 

In the present model the scalar potential is written as 

V 3 3i = V F + V D + V soit , 

where 

V F = -C F = Y J FlF m 



E 

ijk 



+ 



Px / , /a 



+ 



(80) 



+ 



// 2 p/ 2 w 

2 / / . A*X i 72 / / , Mp ,72 / / 

■ylit-j^xi + -f Xt + Y eijkVjPk 2 P Y eijkXjVk 



(81) 



and 
Vb 



-£ D = 1(JW + DI>) = ^(-^ + f/V - X f X + xV + 2pV " WO 2 
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+ 



T^Wj^ - WiKjVj + X^X; - X?A#X$ + pWiPj + PrKjP'j?, (82) 
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finally, 

T/ r 2 t i 2 t i 2 t i 2 /+ / 

Vsoft = -Asoft = "V? 1 T) + m p p ] p + m^X+m^ rj ] rj 

+ m2 P 'P ] p' + - tijkihPiXjVk + kip'iXjrfk 

+ H.c), (83) 

where we have used the scalar multiplets in Eqs. ([57]) and (f69|) . With Eqs. ([HT]) - ([83]) 
we can work out the mass spectra of the scalar and pseudoscalar fields by making 
the usual shift X° — > -^(vx + Hx + iFx)- The analysis is similar to that of Ref. [22] 
and we will not write the constraints equation, etc. By using as input the following 
values for the parameters: sin 2 9w = 0.2314, g = 0.6532, g' = 1.1466; fa = 2, 
f 2 = 10~ 3 ; ki = fcj = 10 GeV; p v = p p = p x = -10 3 GeV; m rj = 15 GeV, m p = 10 
GeV. m p = 244.99 GeV; m X2 = = 10 3 GeV and m p ' = 13 GeV, we obtain the 
masses 

mi « 1702, m 2 « 1449, m 3 « 387, 

m 4 380, m 5 « 361, m 6 ps 130, (84) 

for the scalar sector (all masses are in GeV). Note that the lightest neutral scalar is 
heavier than the lower limit of the Higgs scalar of the standard model, i.e., m# ~ 114 
GeV. For the pseudoscalar sector we obtain 

Mi w 1702, M 2 » 1449, M 3 » 363, 

M 4 w 5, M 5 = 0, M 6 = 0, (85) 

only the two massless pseudoscalars are exact values, i.e., there are two Goldstone 
bosons as it should be. Notice that there is a light pseudoscalar. A carefully study 
shows that there is a Higgs boson satisfies the conditions for SIDM [43]. 

4 Conclusions 

In conclusion, it is a remarkable fact that the 3-3-1 model has an option for self- 
interacting dark matter without the need of imposing any new symmetry to stabilize 
it. We have shown that the Spergel-Steinhardt bound for self-interacting dark mat- 
ter [30] can be realized in the 3-3-1 model with a reasonable choice of the values 
of the parameters. The 3-3-1 model with RH neutrinos provides two Higgs bosons: 
one is scalar or CP-even and another is pseudoscalar or CP- odd particle having 
properties of candidates for dark matter. In difference with the previous candidate 
which introduced by hand, our self-interacting dark matter arises without impose 
new properties to satisfy all the criteria. From the conditions for SIDM one get the 
bounds for scalar Higgs bosons: = 7.75 MeV in the 3-3-1 model with exotic lep- 
tons [16] and 4.7 MeV < nih < 23 MeV in the 3-3-1 version with RH neutrinos [17]. 
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Scalar dark matter candidates have been recently investigated in [44]. It means that 
the candidate for SIDM has not to be introduced ad hoc as in other models [45]. 
This feature remains valid in the supersymmetric 3-3-1 model with right-handed 
neutrinos. The spin ^ exists in the version with exotic neutral lepton. This would 
turn out to be particularly interesting if observations reveal to be in favor of light 
dark matter. This work was supported in part by National Council for Natural 
Sciences of Vietnam contract No: KT - 04.1.2. 
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